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A three-dimensional weak topological insulator (WTI), being equivalent to stacked layers of two- 
dimensional quantum spin-Hall insulators, accommodates massless Dirac electrons on its side surface. A 
notable feature of WTIs is that surface states typically consist of two Dirac cones in the reciprocal space. 
We study the Landau quantization of Dirac electrons of WTIs in a perpendicular magnetic field. It is 
shown that when the magnetic length Ib is much larger than the interlayer distance a, surface electrons 
are quantized into Landau levels according to the ordinary quantization rule for Dirac electrons. It is also 
shown that, with decreasing Ib toward a, each Landau level and its spin state become modulated in a 
nontrivial manner. We demonstrate that this is attributed to the mixing of two Dirac cones induced by the 
discreteness of the layered structure. 


1. Introduction 

Three-dimensional (3D) weak topological insulators 
(WTIs) can be regarded as stacked layers of two- 
dimensional (2D) quantum spin-Hall (QSH) insula¬ 
tors. The stacking direction is designated by the 
weak vector iz = ( 1 ^ 11 ^ 21 ^ 3 ) with i>i, 1 / 2 , and 1 Z 3 being 
weak indices. As a 2D QSH insulator possesses gap¬ 
less excitations only at its edge in the form of a one¬ 
dimensional (ID) helical channel,a WTI accommo¬ 
dates low-energy electrons arising from helical edge chan¬ 
nels only on its side surface. We refer to such surface 
electrons as Dirac electrons since they obey the massless 
Dirac equation in the low-energy limit. Notably, low- 
energy surface states of WTIs typically consist of two 
Dirac cones in the reciprocal space. Although disorder- 
induced scattering between two Dirac cones is not for¬ 
bidden, it does not necessarily extinguish the topological 
nature of WTIs.®“^®^ Several materials have been pro¬ 
posed as possible WTIs.^^“^^) 

A characteristic feature of Dirac electrons shows up 
when a perpendicular magnetic field is applied. Dirac 
electrons are quantized into Landau levels following the 
unique quantization rule that the energy of the nth Lan¬ 
dau level is proportional to This has been exper¬ 

imentally observed in strong topological insulators,^^’^^^ 
as in graphene.Several authors have theoretically 
studied how the Landau levels are formed on a surface 
of strong topological insulators.However, an at¬ 
tempt similar to this has been lacking for Dirac electrons 
in WTIs. In this paper, we study the Landau quanti¬ 
zation of Dirac electrons on a side surface of WTIs in 
the presence of a perpendicular magnetic field. To de¬ 
scribe Dirac electrons in a magnetic field, we mainly em¬ 
ploy an effective 2D Hamiltonian composed of coupled 
ID chains of helical edge channels.^^’^®’^*’^ It has been 


demonstrated^°^ that this effective 2D model can be de¬ 
rived from a 3D Wilson-Dirac Hamiltonian for bulk topo¬ 
logical insulators.We partly employ this 3D Hamilto¬ 
nian to confirm the validity of the 2D Hamiltonian. The 
central question of this study is: does the Landau quan¬ 
tization in WTIs reveal a unique behavior arising from 
the characteristic features of the system under consider¬ 
ation? The answer is yes. If the side surface is infinitely 
long as we assume below, Landau levels are specified by 
the index n and the wave number ky in the longitudi¬ 
nal direction. In the weak magnetic field regime where 
the magnetic length Ib is much larger than the inter¬ 
layer distance a, the behavior of the Landau level is ba¬ 
sically explained within the ordinary quantization rule 
for Dirac electrons. However, a nontrivial behavior ap¬ 
pears with decreasing Is toward a. We find that, when 
Ib/o. becomes sufficiently small, each Landau level and 
its spin state are modulated in an oscillatory manner as 
a function of with period 1 and the modulation 

becomes pronounced with increasing index n. We demon¬ 
strate that this nontrivial behavior is attributed to the 
mixing of two Dirac cones induced by the discreteness of 
the layered structure of WTIs. 

In the next section, we present the effective 2D Hamil¬ 
tonian as well as the 3D Wilson-Dirac Hamiltonian, and 
briefly explain the relationship between them. In Sec. 3, 
we present the massless Dirac equation that is derived 
from the 2D Hamiltonian in the continuum limit. The 
quantization rule for Dirac electrons is derived on its ba¬ 
sis. In Sec. 4, we numerically determine the band struc¬ 
ture of Dirac electrons and calculate the spin expectation 
value of lower quantized levels for several magnetic field 
strengths. We find that each Landau level and its spin 
state become modulated in a nontrivial manner with de¬ 
creasing Ib toward a. In Sec. 5, we present an explana- 
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Fig. 1. Structure of a WTI sample considered in the text; it 
consists of M layers in the z-direction with M > j > 1 and N sites 
in the ^-direction with N > i > 1, while it is infinitely long in the 
y-direction. A magnetic field B is applied perpendicularly to the 
side surface on the j/ 2 -plane. 


tion for the nontrivial behavior. Section 6 is devoted to a 
summary and discussion. We set h = I throughout this 
paper. 


2. Model 


Let us consider a WTI sample that consists of M layers 
in the 2 -direction and N sites in the x-direction, while 
it is infinitely long in the y-direction (see Fig. 1). Let j 
be the index used to specify the M layers, each of which 
is regarded as a 2D QSH insulator. The N sites in the 
a;-direction are specified by the index i. We assume that 
low-energy Dirac electrons emerge on the side surfaces 
of height M in the j/z-plane at i = 1 and If is 
chosen to be sufficiently large, the low-energy states on 
one surface and those on the other surface are mutually 
independent and exactly degenerate. Thus, we hereafter 
consider only the side surface at z = 1. Let |j)^ and |j)^ 
respectively be the basis vectors for right-going and left¬ 
going branches of the ID helical channel arising from the 
jth QSH layer. Under the translational invariance in the 
y-direction, surface electron states are characterized by 
the wave number ky. The effective 2D Hamiltonian for 
Dirac electrons on the side surface is given 


M 

H2D = E \J) 


i=i 


vky 0 

0 —vky 


{j\ 


M-1 







01 + h.c. L (1) 


where \j) = {Ij)-!-, |j) 4 .} and a is the lattice constant, 
which is simply the interlayer distance. In this model, 
the Dirac points appear at {ky,kz) = (0,0) and (0,7r/a). 
For convenience, we define the 2 -coordinate of the jth 
chain as 


M + 1 


( 2 ) 


by setting the origin (i.e., 2 = 0) at the center of the 
system. A perpendicular magnetic field B = {B, 0, 0) is 
introduced in terms of the vector potential A = (0, Ay,0) 


with 


Ay = —Bz. 


( 3 ) 


With this choice of A, the wave number ky remains a 
good quantum number since the translational invariance 
in the y-direction is preserved. The effect of B can be 
taken into account by replacing ky in the Hamiltonian 
with ky + eAy. The Zeeman term is ignored in the main 
part of our analysis, and its effect is briefly discussed in 
Sec. 6. 

To supplement H 2 D, we also use the 3D Wilson-Dirac 
Hamiltonian for topological insulators:^®^ 


H3D 


Mk 

A_Lfc^a 

0 

A||fc_a 

A±kza 

-Mk 

A||fc_a 

0 

0 

A||fc+a 

Mk 

—A±kza 

A||fc+a 

0 

—A±kza 

-Mk 


where k± = k^ ^ iky and 

Mk =mQ + rn2\\(k1 + ky)a^ A m2i.kla^. 


( 4 ) 

( 5 ) 


We discretize the x- and 2 -coordinates and implement 
iLsD on the square lattice in the a; 2 -plane leaving the 
y-coordinate unchanged.In accordance with our as¬ 
sumption that the system is equivalent to 2D QSH insu¬ 
lators stacked in the 2 -direction, we focus on the weak 
topological phase with u = (001). After the discretiza¬ 
tion, the Wilson mass term Mk is modified to 


= mo -I- m2|| {2[1 - cos(fca;a)] + (fcj/o)^} 

+ TO 2 _l 2[1 - cos(fc^a)]. (6) 


The weak topological phase with u = (001) is stabilized 
when the parameters satisfy^^^ 

^211 > ^\rno\ > m2± > ^\mo\- rn2\\, (7) 

where m 2 || > 0 > mg is assumed. This condition fixes 
the sign of the mass term at four symmetric points on 
the kxkz-'p\a,ne with = 0 as follows: 


{ mo <0, k' = (0,0) 

mo -I- 4m2|| >0, fe' = (tt/u, 0) 

mo-I-4m2_L < 0, fe' = (0,7r/a) ’ 

mo -I- 4m2|| -I- 4m2_L >0, fe' = (tt/u, tt/o) 

( 8 ) 


where fe' = {kx,kz). This indicates that, on the side 
surface in the y 2 -plane, the Dirac point appears at 
{ky,kz) = (0,0) and (0,7r/a) in accordance with the 2D 
model. Again, the effect of B can be taken into account 
by the replacement of ky with ky -|- eAy. 

It should be pointed out that the 2D model can be 
derived from the 3D model.The velocities in the 2D 
model are directly related to the parameters of the 3D 
model. 


V = AyO, 


(9) 
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v' = A±a. (10) 

3. Continuum Dirac Theory 

Let us consider the effective 2D Hamiltonian given in 
Eq. (1) in the continuum limit. Note that our model in¬ 
volves the two Dirac cones. Hereafter, the Dirac cones 
centered at {ky,kz) = (0,0) and (0,7r/a) are respectively 
referred to as the first and second Dirac cones. In the 
continuum limit, the effective Hamiltonians and H_ 
describing low-energy states in the first and second Dirac 
cones are given by 

v{ky — eBz) l^v'dz 

±v'dz —i’{ky — eBz) 


H± = 


( 11 ) 


where the vector potential is explicitly included. It is 
convenient to parameterize the strength of B in terms of 
the magnetic length defined by 

The anisotropy of the system is characterized by 

r=^. (13) 

We determine the eigenstates of H± by using the an¬ 
nihilation and creation operators satisfying the commu¬ 
tation relation of [a, al] = 1, defined by 


a = 


al = 


Ib 

V2 

Ib 

71 


^{z - Zc{ky)) + dz 

li 


^{z - Zc{ky)) - dz 


(14) 


(15) 


where Ib = ^/tIb and Zc{ky) = l%ky. With these opera¬ 
tors, H± is rewritten as 


H± = 


\/rv 

V^li 


— (a -I- al) =F (a — al) 
± (a — al) (a -I- al) 


(16) 


The eigenvalues are quantized into the Landau levels as 

'Jrv 


E±n = ± - 


h 


/ / 

-v^ 


(17) 


with n > 0. The wave functions are expressed in terms 
of the eigenstates of the number operator a^a satisfying 


a'^aipn{z) = n-ipniz), 


(18) 


where the explicit forms of the lower three functions are 
given by 


i’oiz) = 

= 


1 


(Z-Zclky))-^ 


\f2 ( z - Zciky) 


_(z-zc(kj^ 

e ^ 


(19) 

( 20 ) 


( 21 ) 

Note that ipniz) depends on ky since it is centered at 
Zc{ky). The eigenstate of with the energy E±n is 
expressed as 

for n > 0 and 

The eigenstate of H- with the energy E±n is expressed 
as 

^ 

for n > 0 and 

The spin state of cj)^ with m = 0, ±1, ±2, • • • is deter¬ 
mined by the expectation value of the Pauli spin matri¬ 
ces, 


(<)'m W\ <)>m) = / ^ 

(26) 

We easily obtain 


(</)±|^|</.±) = (0,0,0) 

(27) 

for m 0 and 


(<(>0 = (±1:0,0). 

(28) 


The above argument provides the quantization rule 
for Dirac electrons in WTIs. Firstly, the energy of the 
nth Landau level is quantized into a value proportional 
to 7^, as in monolayer graphene.Secondly, each 
Landau level is doubly degenerate, reflecting the exis¬ 
tence of the two Dirac cones. Thirdly, the Landau levels 
are not spin-polarized except for the 0th Landau level 
which is polarized to the ±a;-direction. This rule 
is reliable when Ib/o. ^ 1. 

4. Numerical Results 

We numerically determine the band structure of Dirac 
electrons in the presence of a perpendicular magnetic 
field on the basis of the effective 2D Hamiltonian i72D- 
We also calculate the expectation values of z and cr to 
observe whether the quantization rule obtained from the 
continuum Dirac theory is applicable even in a strong 
magnetic field. 

Figure 2 shows the band structure of the system with 
r = v'/v = 0.4 and M = 21. We see that several 
dispersion-less flat subbands appear in a restricted region 
of ky. They correspond to Landau levels. Inside the dis- 


’^oiz) 
■00 (^) 


Ipniz) =F tpn-liz) 
'4^niz) ^ tpn — liz) 


(24) 


■00 (^) 
i^oiz) 


(23) 


'Pniz) ^ pn — liz) 
Ipniz) ± Ipn-liz) 
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Fig. 2. Band structures in the cases of (a) /s/a = 4, (b) /s/a = 
3, and (c) l^fa = 2.4 with r = 0.4 and M = 21. 


played energy range of < 0.5, we can identify the 

Landau levels with the indices n = 0, ±1, ±2, ±3, and 
±4 in the weak magnetic field case of Zs/a = 4. In the 
stronger case of Z^/a = 3, the Landau levels with n = ±3 
and ±4 cannot be identihed. Note that our system is re¬ 
stricted to the finite width of (M — l)a/2 > \z\. Thus, 
roughly speaking, the Landau levels with a small n stably 
exist under the condition of {M — 1)12 —Isla > Z^|fcy|/a 
as ipniz) is centered at Zc = l%ky and its spatial range 
is on the order of In the case of M = 21 and 

lB/a = 'I, the condition is simplified to 6 > l%\ky\/a, in 
good agreement with the result shown in Fig. 2. Outside 
of this region, each energy level shows a dispersion as 
a function of ky. This suggests the appearance of chiral 
edge states. We also see that each Landau level is dou- 




Fig. 3. Expectation values in the case of Zs/a = 4 with r = 0.4 
and M = 21: (a) {z)/a and (b) (ctj) in the 0th Landau level (solid 
line) and in the 1st Landau level (dotted line). 


bly degenerate, reflecting the existence of the two Dirac 
cones. A peculiar feature of the Landau levels shown in 
Fig. 2 is that they deviate from the flat dispersion and 
are modulated in an oscillatory manner as a function of 
Igky/a with a period of 1. This becomes pronounced with 
decreasing Ib and/or increasing the index n. The oscilla¬ 
tory behavior cannot be explained within the continuum 
Dirac theory. 

Figure 3 shows the expectation values of z and CTz in 
the 0th and 1st Landau levels in the case of Ib/o, = 4 
with r = 0.4 and M = 21. We see that {z)/a is propor¬ 
tional to ky in the central region of 6 > l'^\ky\/a. Outside 
of this region, {z)/a approaches the edge position (i.e., 
±10) with increasing \ky\, indicating the appearance of 
chiral edge states. This is also confirmed by the result of 
{az) in the region of l%\ky\/a > 10, indicating that the 
spin is polarized to the ±z-direction (i.e., (cr^) ~ ±1) 
with increasing \ky\. As an edge state going to the pos¬ 
itive (negative) direction is mainly composed of spin-up 
(spin-down) channels near the corresponding edge of the 
system, its spin should be polarized to the z-direction 
(—z-direction). In the central region of 6 > lg\ky\/a, 
{(Jz) vanishes in accordance with the continuum Dirac 
theory. The expectation value of as well as that of ay 
vanishes not only in the 1st Landau level but also in the 
0th Landau level regardless of ky (data not shown). This 
is in contrast to the prediction of the continuum Dirac 
theory (i.e., the 0th Landau level is spin-polarized to the 
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Fig. 4. Expectation values in the case of Zg/a = 2.4 with r = 0.4 
and M = 21: (a) (z)/a and (b) (cr^) in the 0th Landau level (solid 
line) and in the 1st Landau level (dotted line). 


±a;-direction), implying that the corresponding states in 
the first and second Dirac cones are hybridized to form 
spin-unpolarized states. 

Figure 4 shows the expectation values of 2 and az in 
the 0th and 1st Landau levels in the case of Is/a = 2.4 
with r = 0.4 and M = 21. We see that (az) oscillates 
as a function of Igky/a with a period of 1. In the 1st 
Landau level, {z)/a also shows an oscillatory behavior. 
These features cannot be explained within the continuum 
Dirac theory. 

One may think that the 2D model becomes unreliable 
with decreasing Ib and the oscillatory behavior shown 
above is an artifact induced by an erroneous applica¬ 
tion of the 2D model. In the remainder of this section, 
we show that this hypothesis can be ruled out. To show 
this, we calculate the band structure in the case of much 
stronger magnetic fields by using both the 2D model and 
the 3D model. In determining the band structure on the 
basis of iLsD, we consider an infinitely long rectangu¬ 
lar prism-shaped system of height M and width N (see 
Fig. I). Setting M = 21 and N = 20, we numerically 
determine the band structure with the following param¬ 
eters: A_\_/A\\ = 0.4, mo/A|| = —0.6, m 2 \\/A\\ = 0.5, 
”t- 2 _l/^|| = —0.1. Figure 5 shows the band structures 
obtained from the 2D and 3D models in the cases of 
Ib/cl = 2 and 1. The results are almost identical with 
each other. This supports the reliability of H 2 D even in 
a strong magnetic field. 


(a) 




Fig. 5. Band structures in the cases of (a) Zg/a = 2 and (b) 
Zg/a = 1 with r = 0.4 and M = 21, where solid and dotted 
lines, respectively representing the results of the 2D and 3D models, 
almost completely overlap with each other. 


5. Beyond the Continuum Dirac Theory 

In the previous section, we observed the nontrivial 
behavior of the Landau level that cannot be explained 
within the continuum Dirac theory. We can consider that 
this reflects the bare character of ID helical channels con¬ 
stituting surface states of WTIs, as explained below. 

As a transparent example, let us focus on the subband 
corresponding to the 0th Landau level. If the mutual cou¬ 
pling of ID helical channels is ignored, the energy of the 
jth channel disconnected from neighboring ones is 

Cjo- = (^y ~ 

under the choice of A given in Eq. (3), where cr = -I- 
for the spin-up state and a = — for the spin-down state. 
Hence, at 


the energy of the jth channel is ej„ = 0 regardless of 
its spin state while that of the nearest-neighbor chan¬ 
nels is Cjiicr = ^CFva/Pg. Thus, we see that the coupling 
of the jth helical channel to the j -I- Ith and j — Ith 
channels becomes weak with decreasing /s/a. This indi¬ 
cates that, in the vicinity of l^ky/a = Zj/a, the subband 
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state is governed by the jth channel if Zb/o is sufficiently 
small, resulting in nearly linear energy dispersions. This 
behavior is most clearly observed in the case of /b /a = 1 
shown in Fig. 5(b). Generally speaking, the coupling of 
ID helical channels is effectively weakened by an applied 
magnetic field and, as a result, the character of nearly 
disconnected ID helical channels manifests itself when 
Ib/cl is sufficiently small. 

The explanation given above is based on an approach 
from the strong magnetic field limit. An alternative ap¬ 
proach from the weak magnetic held limit is also possible. 
We show below that the nontrivial behavior observed in 
the previous section is explained within the continuum 
Dirac theory by adding a correction that induces the mix¬ 
ing of the two Dirac cones. Our attention is focused on 
surface states located near the center of the system, and 
thus chiral edge states are beyond our consideration be¬ 
low. 

We start with the approximate solutions of the eigen¬ 
value equation for i? 2 Dj dehned by 

\<l>±n) = VaYl ( 31 ) 

j 

kin) = (32) 

j 

where the factor is attached to \4’±n) since it be¬ 
longs to the second Dirac cone centered at {ky,kz) = 
(0, tt/o). In the limit of Ib ^ a, they satisfy H2D\4’±n) = 
Ain kin) as well as the orthonormalization condition of 

(<))±nkin') = ^a,a'^n,n' with a, a' = + Or —. However, 
with decreasing Ib toward a, these equations are satisfied 
only in an approximate sense. Note that the term repre¬ 
senting electron transfer between neighboring chains in 
H 2 B is approximated by the derivative with respect to z 
in the Dirac theory. This approximation becomes worse 
with decreasing Ib ■ Let us consider a correction for this 
approximation. 

For definiteness, we consider the nth Landau level with 
En> 0 consisting of |(/)+) in the first Dirac cone and \(j)~) 
in the second Dirac cone, and evaluate the energy devia¬ 
tion from En by taking account of only the hybridization 
between them. Since (pni^) is the solution of the Dirac 
equation, satisfies 

{H2D - A„) \cj)^) = bVyLf 

j 

X [2adz(l)^{zj) - 4>t{zj+i) + <f)i{zj-i)] , (33) 

where = 1 and P~ = The right-hand side of 
Eq. (33), including the difference between the derivative 
and the finite difference, represents the error of the con¬ 
tinuum approximation. If ) = Sa^a' is assumed, 

we obtain 

(<^±|iL 2 Dkn)=A„ + <5„, (34) 



Fig. 6. ky dependence of the energy of the nth Landau level in 
the case of Ib/o^ = 2.4 for (a) n = 0 and (b) n = 1, where solid 
lines represent the result obtained from iL 2 D while dotted lines 
represent Eq. (41). 



{(j)t\H2r,\4>t) = In 

(35) 

with 



Sn = 

^ \ — l Ty^n—l,n\ 

j’ 

3 

(36) 

In = 

j 

(37) 

for n > 0 and 

s« = 0, 

(38) 



(39) 


j 


where 

/C’™ = IpniZj) [2adz1pm{Zj) - Ipmizj+l) + Ipm^Zj-i)] . 

(40) 

Sn and 7 „ appear to represent the correction arising from 
the discreteness of the layered structure. Note that they 
depend on ky as ipn{zj) is centered at Zc = I'sky. Diago¬ 
nalizing the 2x2 matrix composed of Eqs. (34) and (35), 
we find that the energy is shifted to 

E^{ky) = E^ + Sn{ky)±'y,i{ky). (41) 

In Eig. 6, Eq. (41) is compared with the energy of the nth 
Landau level obtained from H 2 B in the case oUb/o, = 2.4 
for n = 0 and 1. We see that the modulation of the en¬ 
ergy is roughly explained within this treatment. Since 
jn as well as 6n is proportional to the difference between 
the derivative and the finite difference, it becomes large 
with increasing n. This accounts for the observation that 
the modulation becomes pronounced with increasing in¬ 
dex n of the Landau level. The absence of quantitative 
agreement implies that the coupling between the Landau 
levels with different indices is also important. 
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The eigenstate corresponding to is obtained as 

\^t) = VaXl (^i) ± {-^Y<f>n{zj)] , (42) 


which results in 

S+N _ ^ 




2y/2 




(1 + (- 1 )^) iYn{Zj) - 'Ipn-liZj)) 

(1 - (- 1 )^) {i>niZj) + ^n-l{Zj)) \ ’ 


|4>-) 





(43) 


„ r (1 - (-1)^) (YniZj) - 'tpn-liZj)) 
[ (^ + i-^y) (YniZj) + i^n-liZj)) 
for n > 0 and 


II 

+ 0 

\ (i + (-iV) 
[ (i-(-iv: 

1-1 

0 0 

i 

■ (i-(-i)^; 

. ( 1 +(-!)'; 

1 - 1 

0 0 


(44) 


(45) 

(46) 


The above expressions indicate that in |4>+), the spin al¬ 
ternately polarizes to the z-direction in the chain with 
an even j and to the — z-direction in the chain with 
an odd and the direction of spin polarization is com¬ 
pletely reversed in |4>“). Owing to this oscillatory be¬ 
havior of the spin direction, the spin state is not po¬ 
larized owing to self-averaging unless Ib is very small. 
This accounts for the absence of spin polarization in the 
case of Ib/ci = 4 [see Fig. 3(b)]. If Ib decreases to a 
smaller value, self-averaging becomes ineffective and the 
spin tends to polarize to the +z- or — z-direction depend¬ 
ing on whether the corresponding wave function is dom¬ 
inant on the chains with an even j or those with an odd 
j. This accounts for the appearance of oscillatory spin 
polarization in the case oi Ib/ a = 2.4 [see Fig. 4(b)]. 

The alternating behavior of |4>^) is very natural be¬ 
cause iJ 2 D allows a spin-up (spin-down) ID channel to 
couple with only neighboring spin-down (spin-up) chan¬ 
nels. That is, each eigenstate of H 2 D should be repre¬ 
sented by the superposition of helical channels with their 
spin direction alternately flipping from chain to chain. 
This is precisely the feature represented by Eqs. (43)- 
(46). 

Let us consider the strong magnetic field regime where 
Is /a is sufficiently small, focusing on the 0th Landau 
level. The spatial range of the superposition is on the 
order of Ib, so |$^) is dominated by the channel nearest 
to the center of the wave function. This indicates that 
our argument based on the continuum Dirac theory is 
continuously connected to the picture of nearly discon¬ 
nected helical channels argued in the beginning of this 
section. 


6. Summary and Discussion 

In this paper, we have studied how Dirac electrons on a 
surface of weak topological insulators are quantized into 
Landau levels in the presence of a perpendicular mag¬ 
netic field B. When B is sufficiently weak with Ib/o- ^ 1, 
the ordinary Landau level structure with dispersion-less 
flat subbands is observed in accordance with the quan¬ 
tization rule for Dirac electrons. The spin state of each 
level is not polarized except for chiral edge states ap¬ 
pearing near the edges. With increasing B, duS Ib/ci ^ 1, 
each flat subband becomes oscillating as a function of 
l%ky/a with period 1. The corresponding spin state also 
becomes oscillating with the same period. We have shown 
that this nontrivial behavior is attributed to the mixing 
of the two Dirac cones induced by the discreteness of the 
layered structure. Alternatively, it can be regarded as 
a manifestation of the character of a ID helical channel 
nearly disconnected from neighboring ones by a magnetic 
field. We have also shown that these two explanations are 
consistent with each other. 

Here, we discuss the possibility of detecting the non¬ 
trivial behavior of the Landau level in actual experimen¬ 
tal situations. As its simplest consequence, we expect 
the broadening of each Landau level, which will be de¬ 
tectable if its width is sufficiently large. When a mag¬ 
netic field of B [T] is applied, the magnetic length is ex¬ 
pressed as Zb = 25.65/'/B [nm]. While the magnetic field 
of H = 40 T yields Ib ~ 4.06 nm, the interlayer distance 
of the possible weak topological insulator KHgSb^^^ is 
a ~ 1.9 nm. These parameters yield Ib/cl k, 2.1. From 
Fig. 2(c) for the case of Ib/cl = 2.4 and Fig. 5(a) for 
the case oi Ib/cl = 2, we see that this value would re¬ 
sult in sufficiently large level broadening for its detec¬ 
tion. Similar broadening would also appear in monolayer 
graphene.However, an extremely strong magnetic 
field is necessary to detect it as the lattice constant of 
graphene is on the order of 1 A. 

In the presence of disorder, a further consideration 
is needed as disorder also induces level broadening. 
Roughly speaking, the broadening due to disorder is pro¬ 
portional to oc ^/~B and is independent of the index 
Hence, the broadening due to disorder slowly in¬ 
creases with increasing B, while the oscillation-induced 
broadening rapidly increases in the strong magnetic field 
regime, as can be seen from Figs. 2 and 5. Furthermore, 
the former is independent of the index n while the latter 
becomes large with increasing n. We expect that these 
qualitative differences enable us to detect the oscillation- 
induced broadening as long as the disorder is not so 
strong. 

So far, the Zeeman effect on Landau levels has been ig¬ 
nored in our argument. Here, we briefly consider it within 
the continuum Dirac theory. If we take into account the 
Zeeman term for a magnetic field B = (H,0,0), the 
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Hamiltonian H± becomes 


H± 


v{ky — eBz) + A 

±v'dz + A —v{ky — eBz) 


with 


A = 


1 


(47) 


(48) 


where /ie and 9 are respectively the Bohr magneton and 
the effective g factor. In the presence of the Zeeman term, 
the energy of the nth Landau level is modified as 


where 



(49) 


(50) 


If u = 3 eVA is assumed, we hnd that A r; O.Ollg/v^ at 
i? = 40 T. The above argument implies that the Zeeman 
effect induces only a small shift of each quantized level 
and does not strongly alter the qualitative behavior of 
the magnetic field effect on Landau levels. 

Let us finally examine the stability of surface states 
against an increase in the magnetic field. For definite¬ 
ness, we consider the 0th Landau level in the vicinity 
of llky/a = Zyla. The corresponding wave function is 
centered at the jth chain and its spatial range is on the 
order of Z_b. A plausible criterion to ensure its stability is 
that all helical channels within the distance of Ib from 
the jth chain stably exist. Equations (6) and (8) indicate 
that this is ensured under the condition of 

rho + m2\\ <0 (51) 


with mg = max{TOo, Wo + 4 to2_l}, where m 2 || > 0 > mg 
is assumed. It is convenient to define the characteristic 
length as 1%/a = (m 2 ||/|mo|)^/^. If Is > 1%: condi¬ 
tion (51) is satisfied, resulting in the stabilization of the 
surface state. However, this does not mean that the sur¬ 
face state completely disappears when 1% > Ib- Even 
in this regime, the helical channel in the jth chain per¬ 
sists and its energy is zero at l%ky/a = Zj/a^ but helical 
channels in its neighboring chains within the distance of 
Ib are partly or fully destabilized and are replaced with 
bulk states. As a consequence, the persistent jth heli¬ 
cal channel is inevitably coupled with such bulk states. 
Hence, although the resulting state contains the jth he¬ 
lical channel, it cannot be regarded as a pure surface 
state. 
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